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Weak lensing predictions for modified gravities at 
non-linear scales 
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ABSTRACT 

We present a set of predictions for weak lensing correlation functions in the context of 
modified gravity models, including a prescription for the impact of the nonlinear power 
spectrum regime in these models. We consider the DGP and f(R) models, together 
with dark energy models with the same expansion history. We use the requirement that 
gravity is close to GR on small scales to estimate the non-linear power for these models. 
We then calculate weak lensing statistics, showing their behaviour as a function of scale 
and redshift, and present predictions for measurement accuracy with future lensing 
surveys, taking into account cosmic variance and galaxy shape noise. We demonstrate 
the improved discriminatory power of weak lensing for testing modified gravities once 
the nonlinear power spectrum contribution has been included. We also examine the 
ability of future lensing surveys to constrain a parameterisation of the non- linear power 
spectrum, including sensitivity to the growth factor 7. 

Key words: Gravitation; Gravitational Lensing; Cosmology: Theory 



1 INTRODUCTION 

Consistent observational evidence from various cosmological 
probes shows that the Universe is currently undergoing a 
period of accelerated expansion. The observed expansion 
history can be explained using some form of dark energy 
or a cosmological constant; however this cosmological 
constant cannot be explained with current particle physics 
due to its very small value. An alternative approach is to 
invoke a modification of gravity; there are many different 
ways that gravity and/or the equation of state of the 
dark energy can be modified to allow for the expansion 
history observed. This makes it impossible to differen- 
tiate between the effects of modified gravity and dark 
energy by measuring the background expansion history 
alone. However, modifying gravity also produces a distinct 
growth rate of structure; thus the expansion history and 
growth history together can be used to distinguish be- 
tween various models of gravity. This consistency relation 
to test GR has been propo s ed an d explored by ma ny 



Pan-STARR^ and LSSlB and future space surveys such 
as Euclic0, will allow a combination of growth of structure 
and expansion history to be probed to considerably higher 
precision, which will allow many gravity models to be 
excluded. 

There has been a great deal of work showing how to 
use weak lensing to discriminate between different gravity 
models; however this has been restricted to probing the lin- 



ear r e gime of the ma t ter power spe c trum (Afshord i et al 
20081 ; ISchmldl 120081; feong fc Dorel 120081; iThomas et al 



20081 ; jTsuiikawa & Tatckawa 2008; IZhao et al.l l2009al lbh or 



to test IjK nas been propo s ed an d explored by many 
papers (|Uzan fc Bernardeaul [200ll; iLue et alj [2004al|3; 
Ishak et al . 2006; K unz fc Saponejl2007l; Ichiba fc Takahashil 



20071; IWang et al.1 I2007T: 



Jain fc Zhand 120081; 



Bertschinger fc Zukinl 2008; 
Daniel et al.l 120081 ; ISong fc Kovamal 



2009). Upcoming weak lensing surveys such as DE^3, 



uses methods tha t do not obta i n GR at small scales 
ilKnox et all 120061; iHeavens et all 120071 ; lYamamoto et al.1 
120071 ; lAmendola et al ]|2008l ). The non-linear regime provides 
much of the power for lensing and can be most easily probed 
by current and upcoming lensing surveys. This paper exam- 
ines the effect of including the non-linear regime in modi- 
fied gravity lensing predictions, including the small-scale GR 
limit, to see how useful weak lensing will be overall when try- 
ing to determine the correct model of gravity. First we look 
at DGP and f{R) gravity models as examples, and inves- 
tigate weak lensing's ability to differentiate between these 
models and dark energy models. We then take a more phe- 
nomenological point of view, by parameterising the shape of 
the matter power spectrum and examining the sensitivity of 
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weak lensing observables to changes to the matter distribu- 
tion when the expansion history is the same for each model 
considered. Using these parameters we show how strongly 
a ground-based survey similar to DES and a space-based 
survey such as Euclid will be able to discriminate between 
different growth histories with identical expansion histories. 

This paper is organised as follows. In section[2]we briefly 
describe the DGP and f(R) models of gravity and how they 
compare with dark energy models. We describe how we cal- 
culate matter power spectra for these models, including the 
GR small-scale limit. We also describe how we proceed to 
calculate weak lensing observables from these power spectra. 
In section[3]we present the resulting lensing correlation func- 
tions, including realistic errors for future surveys taking into 
account shape measurement noise and cosmic covariance. In 
section [4] we take the alternative approach of parameterising 
the non-linear power spectrum, and we investigate how sen- 
sitive weak lensing is to these parameters which go beyond 
the usual growth parameter. We present our conclusions in 
section [5] 

Throughout this paper we will use a flat cosmology with 
the WMAP5+SNe+BAO best fit cosmological parameters, 
which are determined by the background evolution of the 
Universe. We use a ACDM background for both ACDM and 
f(R), in which case we take n s = 0.96, h = 0.71 , Q m = 
0.27±0.02 and ct 8 = 0.81±0.03 (|Komatsu et alj2009h . When 
we use a DGP b ackground, we ha ve n s = 0.998, h — 0.66, 
fi m = 0.26±0.02 (|Fang et alJl200Sh giving a a s = 0.66±0.03 
for an equivalent ACDM model. 



2 LENSING IN DGP AND f(R) MODELS 

2.1 Modified gravity power spectra 

As we have already mentioned, there are two key phenomena 
to model in any gravity in order to calculate the matter 
power spectrum: the expansion history, quantified by the 
evolution of the Hubble parameter, and the growth history, 
quantified by the evolution of density perturbations S in the 
Universe. 

For ACDM, the expansion history is given by the Fried- 
mann equation 

H 2 = ^ + n A , (i) 

where H = ^/aH , a is the scale factor and Ho is the 
present day Hubble constant. 

The growth history is described by the density per- 
turbation evolution equation together with the Friedmann 
equation. At this point we will limit ourselves to the regime 
where density perturbations evolve linearly. In this regime 
we have 



3 H' \ , _ 3G ff fim „ 
a H ) ~ 2H 2 a 2 a 3 ' 



(2) 



where primes denote differentiation with respect to a. This 
equation is valid for both dark energy and modified grav- 
ity models, where G c ff is the effective gravitational constant 
normalised by the gravitational constant G; hence G e g = 1 
for dark energy models, while for modified gravity models 

G cff = l + ^, (3) 



where (3 is determined by the model. 

In this paper, we consider DGP (|Dvali et al.ll200oh and 
f(R) as examples of modified gravity models, as the non- 
linear power spectra have been studied in great detail in 
these two models using perturbation theory and N-body 
sim ulations. For some review s of modified gravity model s 



Kov 



Noiiri fc Odintsovl l|200rj) ; [5 urrer fc Maartena ([2008) ; 
' l|2008l V~ 



In DGP, spacetime has five dimensions, while we live 
on a 4D brane in the 5D bulk. Standard Model particles 
are bound on the 4D brane, as is gravity on small scales; 
however on large scales gravity leaks off the brane causing 
late time acceleration. The scale of the transition from 4D 
to 5D gravity is governed by the crossover scale, r c = (1 — 
fim) -1 . The extra dimension contributes a further term to 
the Friedmann equation whose amplitude is governed by r c : 



The growth history is also altered, 
ama fc Maartensll2006l ) 



2Hr c 1 + 



aH' 
~3H 



giving 



(5) 



In f(R) gravity models the Einstein-Hilbert action is mod- 
ified to include an arbitrary function of the Ricci scalar, 
R. In this study we use an f(R) function of the form 
|Hu fc Sawickill2007bl ) 



/ = -6«a-^/r , 



(6) 



where R is the Ricci scalar, Ro is the present day Ricci scalar 
and /r = ^n\ R _ R ■ We use |/r | = 10 -4 , which has been 
found to fit with cluster constraints l|Schmidt et al . 200Sj),to 
give a background evolution which is approximately ACDM 
to sub-percent level. This allows us to use the ACDM Fried- 
mann equation and only alter the density evolution e quation 
|Lue et al.ll2004bl ; [Zhang 2006]; iKovama et aill2009l ) with 



P = l + 



1 /a 
3c^ Vfc 



(7) 



where k is the dimensionless wavenumber defined as 
k(c/Ho), k is the wavenumber and c is the speed of light. 

For modified gravity to agree with solar system obser- 
vations it must approach a GR solution on small scales. This 
means that the non-linear power spectrum must be an inter- 
polation of the modified gravity non-linear power spectrum 
with no mechanism to obtain the GR result on small scales, 
Pnon-GR.(fc, z), and the GR non-linear power spectrum with 
the same expansion history as the modified gravity model, 
Pgr(^iz )- A fitting formula fo r this interpolation was pro- 
posed bv lHu fc Sawickil (|2007aT ): 



P(k, 



(k,z) +c ril {z)Y?{k,z)P GK {k,z) 
l + c nl (z)^(k,z) 



(8) 



where E (k,z) picks out non-linear scales and c n \(z) deter- 
mines the scale at which the power spectrum approaches the 
GR result as a function of redshift. 

In this paper, we use the fitting fomul ae for Y?(k,z) 



and c n i(z) obtained by perturbation theory ( Koyama et al 



2009) and confirme d by N-body simulations (jOvaizu et al 
200alSchmidll2009h . 
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Figure 1. Matter power spectrum for ACDM, DGP and f(R) at 
z=0. 



Figure 2. Relative difference between matter power spectra for 
ACDM and f(R) at z = for different ai- 



Z 2 (k,z)= (j^P lin (k,z)j , c nl (z) = A(l + z) a2 . (9) 

where pi n (fc,z) is the modified gravity linear power spec- 
trum. The non-linear power spectrum for both t he P n0 n-GR 
and Pgr is found using the ISmith et all (|2003l ) fitting for- 
mula from the linear power spectrum. For DGP, A = 0.3, 
a\ = 1 and ai = 0.16 and for f(R) with /r = 10~ 4 we 
use A = 0.08, ai = 1/3 and a 2 = 1.05 for ^ z < 1. It 
should be noted these values are not valid for all f2 m and as ■ 
However, in DGP, these values depend on f2 m and as very 
weakly, so within our priors for Q m and as we can assume 
the values are constant. 

We should also emphasise that these fits are confirmed 
only up to k = lfe/Mpc due to the lack of resolution in N- 
body simulations, so we are extrapolating the fits beyond 
this regime. Clearly it is necessary to check the validity of 
this extrapo lation using N-body si mulations with higher res- 
olution (see ISchmidt et al.l ((2008) for a different approach 
using the halo model). However, since the modified gravity 
power spectrum should approach the GR non-linear power 
spectrum with the same expansion history, and since the fit- 
ting formula (JS| ensures this, our extrapolation is justified. 

In applying this formalism, we found that although 
f(R) fits the N-body results at small k, it failed to con- 
verge with ACDM at larger k if ot\ = 1/3. This is due to the 
strong scale dependence of the linear power spectrum, such 
that Pnon-GR deviates from Pgr strongly on small scales 
and equation (JHJ with ot\ — 1/3 fails to converge with Pgr. 
Thus, we also consider Qi = 1 and ai = 2 cases for f(R) 
which have more physical behaviour at high k. 

Since we are interested in how sensitive weak lensing is 
to different growth histories with the same expansion his- 
tory, we will also consider a quintessence cold dark matter 
(QCDM) model. In this case, the equation of state of the 
dark energy is altered to match the expansion history of 
DGP, while the density perturbation evolution equations are 
the same as ACDM. 

We show examples of the resulting matter power spectra 
in Figure[T] f(R) models show the scale dependent enhance- 
ment of the power spectrum in the linear regime compared 
with ACDM. For ai = 1/3, which fits N-body results well 
up to k = 1/i/Mpc, the power spectrum fails to converge 




1.0 10.0 
k[h/Mpc] 

Figure 3. Relative difference between the matter power spectra 
for DGP and the QCDM model at z = 0. 



with ACDM. On the other hand, the power spectrum with 
Qi = 2 shows clear convergence; this is shown more explic- 
itly in Figure [2] 

We also show a comparison between DGP and QCDM 
power in Figure [3] including our non-linear prescription. In 
the linear regime the DGP power spectrum receives scale 
independent suppressions, but it converges to the QCDM 
power spectrum on non-linear scales due to our inclusion of 
the GR asymptote. 



2.2 Weak Lensing 

The mass distribution described by the matter power spec- 
trum deflects light all the way along the path from source 
to observer, so any changes in the matter power spec- 
trum alter the observed image distortion of galaxies. Light 
bundles are transformed by a shear with two components, 
7 = 7i + ?72, and an isotropic dila tion, the convergence k 
fe.g. lBartelmann fc Schneiderl200ll ). The Jacobian mapping 
from the unlensed image to the distorted image is given by 



.4 



1 — ft — 7i 

—72 



-72 

1 — K + 71 



(10) 
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Figure 4. Redshift distributions used for survey predictions: for 
ground-based survey with z m = 0.825, and for Euclid with z m = 
0.91. 



Since each galaxy has some unknown intrinsic shape, the 
amount of lensing cannot be estimated by a single source; 
however one can correlate the shear estimators of many 
sources, in which case randomly oriented intrinsic elliptici- 
ties will average out, leaving the gravitational shear signal. 
This will not succeed if galaxy ellipticitie s are physically 
align ed, which they are to some degree (e.g. lOkumura et al.l 
2009); however, for our purposes we will assume that the 
resulting physical correlation signal can be removed, leaving 
only the lensing signal. Here we will concentrate on the con- 
vergence correlation function, equal to the sum of the shear 
correlation functions, which is related to the convergence 
power spectrum by 



C K (0) = 



dl^P K (l)Jo(W), 

Z7T 



(11) 



where 9 is the angular distance between the correlated 
sources and I is the angular wavenumber. The convergence 
pow er spectrum is related to the ma tter power spectrum by 
(e.g. iBartelmann fc Schneiderll200ll ) 



with 



d X W( X ) 



W( X ) = / dx'G(x') 1 



(12) 



(13) 



where \ is comoving distance, xh is the comoving distance 
to the horizon and G(x) is the normalised distribution of 
the sources in comoving distance, corresponding to a red- 
shift distribution. Equation (|13[1 is valid for flat cosmologies, 
which are all that are considered in this paper. 

We will calculate results for realistic notional surveys: a 
ground-based survey similar to that of the Dark Energy Sur- 
vey (DES), and a space-based survey such as that of Euclid, 
using redshift distributions shown in Figure [4] the redshift 
distribution for our ground-based survey w a s cho sen to be 
the same as for CFHTLS given bv lFu et all (|2008l ) giving a 



This can be ext ended to cross-cor r elate sources at dif - 
ferent redshifts (e.g. iBacon et alj|200d iMassev et al.ll2007T) 
to obtain 



XH Pg ( i y) 

d X W 1 (x)W 2 ( X ) V V , (14) 



where Wi include the galaxy distributions Gi appropriate 
for the ith redshift bin. This equation together with equa- 
tion (|11[) relates the matter power spectra from our gravity 
models to the predicted lensing signal; we will now use these 
tools to calculate lensing predictions for our models. 

It should be noted that in our analysis we have not 
included the effects of baryons, which are known to have 
an effect on the matter power spectrum for k ^ lhMpc~ 
as sho w n inlWhitd h20o4): IZhan fc Knoxl d2004 ): Jjing et all 
(|2006h : iRudd et all fcOOgft : iHearin fc Zentneil (|2009l ) For 
non-radiative gas simulations this changes the amplitude of 
the matter power spectrum by a few percent, however if gas 
cooling and star formation are included this effect could be 
considerably larger. Therefore for a full lensing analysis in 
the non-linear regime these effects must be included as well, 
but since we must recover GR on small scales, results from 
gas simulations for GR can be used to refine the models we 
present. 



3 RESULTS 

We calculate the convergence (combined shear) correlation 
function of equation (|lip for all of our models, and esti- 
mate measurement errors due to intrinsic ellipticity, for the 
notional ground-based and Euclid surveys using bins with 
error CT shapc = %/2<7^ / y/Np^teiO, A0), where <r 7 = 0.3. The 
errors were estimated using 13.3 galaxies arc min -2 and a 
survey area of 5000 square degrees for our ground-based sur- 
vey (as is appropriate for DES), while for Euclid we use 
35 galaxies arc min -2 and 20000 square degrees. The co- 
variance matrix for the intrinsic ellipticity nois e is diagonal 
for b ins in redshift and angular separation (c.f. IBacon et al.l 
l2003h . 

We also include the covariance due to sample variance 
due to the cosmic matter distributio n, C CO a, which is est i- 
mated using the Horizon simulation l|Tevssier et al.ll2009l ). 
3-D convergence maps were calculated from the 3-D over- 
density field, for 75 patches of area 2 square degrees; con- 
vergence correlation functions were then measured in each 
patch. The covariance between the resulting patch correla- 
tion functions was measured as an estimator of the true co- 
variance, in 8 angular separation bins logarithmically spaced 
from V to 90' and in 3 redshift bins (leading to 6 redshift 
pair bins). The diagonal elements of the covariance matrix 
for mean correlation functions are measured to be approx- 
imately 10~ n — 10~ 9 per square degree, making the sam- 
ple covariance the dominant source of error for larger an- 
gles and higher redshifts for both our ground-based survey, 



with diagonal element values of 10 



10" 16 and Euclid, 



10" 



10 . These should be compared with shape noise 



covariance contributions of 10 15 — 10 11 for ground-based 
and 10" 16 -10" 12 for Euclid. The covariances were included 
tribution given bv lHawken fc Bridle! l|2009h giving z m = 0.9. in our x 2 estimations using (c.f. lHartlap et ai1l2007l ) 



median redshift, z m , of 0.825 and for Euclid we used the dis- 
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= E (*-*)(- 



n G — 1 



Ccos H~ ^shape 



n b 



where d is the 'data', here the fiducial ACDM correlation 
function in redshift and angular separation bins; t is the 
alternative gravity model correlation function in those bins, 
n — 75 is the number of realisations of correlation functions 
used in the calculation of C CO s and rib = 48 is the total num- 
ber of bins in angular separation and redshift. Note that we 
use the sample covariance estimate from Horizon (which fol- 
lows ACDM) for both ACDM and QCDM cases; the QCDM 
error bars should therefore only be considered as the correct 
order of magnitude. 

We calculate for each of our models the differ- 
ence in x 2 between the modified gravity model and a 
dark energy model (either ACDM or QCDM), applying 
WMAP+SNe+BAO priors. Note that for ACDM and J (R), 
we used the ACDM background dKomatsu et al. J |2009l) and 
for D GP and QCDM, the DGP background (|Fang et all 
120081 ) was used (see §1). 

Figure [5] shows example results for our ground-based 
survey and Euclid using the central cosmological parame- 
ter values for WMAP+SNe+BAO described in §1; this is 
for the 2-D projection case where we have not divided the 
catalogue tomographically. We see from figures (a) and (b) 
that the difference between models is substantially greater 
in the nonlinear regime (9 < 30') than in the linear regime 
(6 > 30'), as is the amplitude of the signal. As (c) and (d) 
show, it is also the case that the linear correlation function 
is small in the low-# regime, if nonlinear corrections are not 
included. 

We present the \ 2 differences between the modified 
gravities and fiducial dark energy models in Table Q] for 
the 2-D (non-tomographic) cases including non-linear power. 
We see that there is indeed strong discriminatory power be- 
tween modified gravity models and ACDM with the notional 
ground-based survey; the precision of Euclid is even more 
impressive. 

We also compare the constraints on DGP and a QCDM 
model of the same expansion history (i.e. a DGP back- 
ground). The correlation functions for these models are 
shown in Figure [()] One can either consider a QCDM model 
with cosmological parameters equal to their central values in 
a fit to WMAP+BAO+SNe, or more realistically the best 
fit QCDM model to the DGP model obtained by varying 
fi m and og- We see that there is a choice of Q. m and as that 
make the QCDM and DGP models virtually indistinguish- 
able. This is confirmed by the bottom row of Table [T] which 
shows that the difference in \ 2 f° r DGP and this QCDM is 
insignificant. This is clearly partly due to the existence of a 
QCDM model with rather similar growth to the DGP, but 
also because of the low amplitude of the DGP correlation 
function, with the result that the error bars are larger in 
proportion to the signal than for other models. 

The power of future surveys to discriminate between 
gravity models is borne out by the tomographic results. Ex- 
amples of these are shown in Figure [7] where we see the 
different redshift evolutions and amplitudes of the signal in 
the different gravities. Table [5] confirms that using the red- 
shift information affords us better discrimination between 
dark energy and modified gravity models in every case, by a 



Fiducial 


Modified 


Ground-based 


Euclid 


Model 


gravity 


Ax 2 


Ax 2 




DGP 


4 x 10 3 


3 x 10 4 


ACDM 


/OR), 0l = 1/3 


500 


6 x 10 3 


/OR), ai = 1 


200 


2 x 10 3 




/(-R), «l = 2 


10 


500 


QCDM 


DGP 


0.5 


3 



Table 1. Ax 2 for DGP and f(R) using errors from our ground- 
based survey and Euclid, with no redshift information, and us- 
ing priors from WMAP+SNe+BAO. The top section shows re- 
sults compared to ACDM, while the bottom row is compared to 
QCDM. 



Fiducial 


Modified 


Ground-based 


Euclid 


Model 


gravity 


Ax 2 


Ax 2 




DGP 


6 x 10 3 


7 x 10 4 


ACDM 


/OR), "i = 1/3 


600 


8 x 10 3 


/OR), ai = 1 


300 


3 x 10 3 




f(R), <*i = 2 


60 


1 x 10 3 


QCDM 


DGP 


0.5 


5 



Table 2. Same as Tabic [T] but using tomographic information. 
In each case we have redshift bins of width Az = 0.4 between 
z = 0.3 and 1.5. 



factor of 50 to 100%. Because of this, we will only consider 
tomographic results from now on in the paper. 

Table [3] shows the impact of including non-linear power 
on our ability to discriminate between modified gravities. 
Comparing these results with Table [2] we can see the im- 
provement that measurements from the non-linear regime of 
the correlation functions provide. The improvement is very 
substantial, amounting to an order of magnitude in \ 2 dif- 
ference. 

It is important to note that using only the lSmith et al.l 



(2003) formula, without the GR asymptote, causes an over- 
estimation in our ability to discriminate between modified 
gravity and dark energy models as shown in Table 2] This 
can amount to up to a 90% difference in A\ 2 for some mod- 
els, due to the difference in power at small scales that is 
present when there is no attempt to recover GR. This shows 
the importance of careful modelling of the nonlinear regime, 
including the appropriate small-scale GR limit. 



4 PARAMETERISATION OF THE POWER 
SPECTRUM 

The sensitivity of lensing to changes in the matter power 
spectrum will be very important in determining the correct 



Fiducial 


Modified 


Ground-based 


Euclid 


Model 


gravity 


Ax 2 


Ax 2 


ACDM 


DGP 

f(R) 


500 

3 


3000 
20 


QCDM 


DGP 


0.2 


2 



Table 3. Ax 2 if only linear power is included for 8 = 30' — 
90', for 0.4 redshift bins between 0.3 and 1.5 using priors from 
WMAP+SNe+BAO. 
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(a) Including non-linear effects for sources with z m = 0.825, with (b) Including non-linear effects for sources with z m = 0.9 with 
ground-based survey errors Euclid errors 
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(c) Not including non-linear effects for sources with z m = 0.825 (d) Not including non-linear effects for sources with z m = 0.9 
with ground-based survey errors with Euclid errors 

Figure 5. Correlation function predicted for ACDM, DGP and f[R) with error estimates for ground-based survey and Euclid. Models 
are for the central cosmological parameter values fitting WMAP+BAO+SNe described in §1, using the ACDM background (for ACDM 
and f(R)) and the DGP background (for DGP). 



Fiducial 
Model 



Modified 
gravity 



Ground-based 
% difference 



Euclid 
% difference 



DGP -3% 

f(R), ai = 1/3 -40% 

f(R), «i = 1 -70% 

f(R), ai = 2 -90% 



-3% 
-40% 
-80% 
-90% 



QCDM 



DGP 



-70% 



-80% 



Table 4. Percentage difference in A\ 2 if the Smith et al. (2003) 
formula is used with no attempt to fit GR at small scales, com- 
pared to using the Hu & Sawicki fitting formula. All results arc 
tomographic with WMAP+SNe+BAO priors as before. 



theory of gravity or dark energy in the near future. In this 
section we will therefore parameterise the non-linear power 
spectrum, in order to more fully understand what aspect 
of the power spectrum it is which lensing surveys will be 

sensitive to. 

We use the growth factor 7 l|Linderll2005l ) as is used in 
lAmendola et al.l 1 20081 ). but we also include the parameters 
used in the Hu and Sawicki fitting formula (Equation [8] and 



[9]). In the formalism of iLinderl (|2005l ) the growth history, 
g(a), is given by 



g(a) = exp 



/ fi m V 



da 
a 



(16) 



where 7 is set by the model. This parameterisation can- 
not model all theories of gravity, since it does not allow for 
growth histories which have k dependency, such as f(R)- 
It is also only valid for gravity models where the combina- 
tion of $ + ^ i s the same as in GR, w hich is true for DG P 
|Kovamall200ll) a nd f (R) fo r f R0 < 1 (|Ovaizu et alj|2008l ). 
Figures 8(a) and 8(b) demonstrate the dependence of 



the parameters on one another when fitting weak lensing 
predictions for varying 7, A, ai and a-z to a ACDM fidu- 
cial model when Q m and as are fixed at the central values 
fitting WMAP+BAO+SNe. The slight widening in the 7 
constraint as A, ai and «2 increase is due to being able to 
recover ACDM at non-linear scales by increasing A and ai 
as 7 varies. This means that the constraint on 7 degrades 
slightly by including the parameters in the Hu and Sawicki 
fitting formula (A, ai and 02). The constraint obtained by 
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1 10 100 1 10 100 

9 [arc min] 9 [arc min] 

(a) Including non-linear effects for sources with z m = 0.825 with (b) Including non-linear effects for sources with z m = 0.9 with 
ground-based errors Euclid errors 

Figure 6. Correlation function predicted for the QCDM model with the expansion history as DGP and DGP with error estimates 
for ground-based survey and Euclid. The solid lines show the correlation function for the QCDM model for the central cosmological 
parameter values fitting WMAP+BAO+SNe, using the DGP background. The dashed line shows the best fit QCDM model to the DGP 
model obtained by varying f2 m and ag. 




(a) ACDM for z = 0.3 - 0.7 redshift bin (b) DGP for z = 0.3 - 0.7 redshift bin (c) f(R) for z = 0.3 - 0.7 redshift bin 




(d) ACDM for z = 0.7 - 1.1 redshift bin (e) DGP for z = 0.7 - 1.1 redshift bin (f) f(R) for z = 0.7 - 1.1 redshift bin 



Figure 7. Correlation function predicted for ACDM, DGP and f(R) with error estimates for ground-based survey at different z using 
redshift bins with width Az = 0.4. 



marginalising over all fl m and as shown in Figures 9(a) and 



9(b) shows that the constraint for 7 for a ACDM fiducial 
model is very good, as shown in Table measuring 7 within 
20% of its value for the ground-based survey and within 5% 
for Euclid, while the other parameters are difficult to con- 
strain. 

A better constraint on the parameters can be found for 
a growth history that is not ACDM, such as DGP, as shown 



in Figures 9(c) and 9(d) This provides a better constraint 
on A, Qi and a.2, but the constraint on 7 is not as tight, as 
shown in Table [6] measuring 7 within 30% of its value for 



the ground-based survey and within 12% for Euclid. This is 
due to the degeneracy between 7 and the other parameters 
in this instance. These degeneracies can be seen more clearly 
before the results are marginalised over Q m and as as shown 



in Figures 8(c) and 8(d) The large dependence on the other 



fitting parameters demonstrates that care should be taken 
when predicting 7 constraints using this parameterisation. 

One might think then that it is better not to include 
non-linear scales and constrain only 7 on linear scales. How- 
ever, there is substantial extra signal coming from the non- 
linear regime. In fact with our parameterisations, Tables [5] 



8 Emma Beynon, David J. Bacon and Kazuya Koyama 



0.2 0.4 0.6 0. 




0.50 0.52 0.54 0.56 0.58 0.60 

y 




0.55 
7 



(a) ACDM fiducial model with ground-based errors 



(b) ACDM fiducial model with Euclid errors 





0.55 0.63 0.71 0.79 0.87 0.95 



0.60 0.65 0.70 0.75 



(c) DGP fiducial model with ground-based errors 



(d) DGP fiducial model with Euclid errors 



Figure 8. Constraints on 7, ai, ct2 and A from our ground-based survey and Euclid, using 0.4 rcdshift bins between 0.3 and 1.5 for 
the central cosmological parameter values fitting WMAP+BAO+SNe described in §1. The light grey contours show the 68% confidence 
limits and the dark grey show the 95% confidence limits. 



Survey 



Our 

paramctcrisation 



Smith et al. 
Linear (2003) 



Ground 68% 
-based 95% 
68% 
95% 



Euclid 



0.10 
0.24 
0.030 
0.069 



0.23 
0.42 
0.12 
0.23 



0.091 
0.18 
0.026 
0.051 



Table 5. The 68% and 95% confidence limits for the growth factor 
7 obtained for our parameterisation with ACDM as the fiducial 
model compared to those obtained using only linear scales and 
compared to the constraint from using Smith ct al. 2003 to model 
the non-linear. These are marginalised over fl m , erg, A, a\ and 
Q2- 



Survey 



Our 

paramctcrisation 



Smith et al. 
Linear (2003) 



Ground 68% 
-based 95% 
68% 
95% 



Euclid 



0.22 
0.59 
0.082 
0.12 



0.38 
0.68 
0.20 
0.39 



0.25 
0.48 
0.052 
0.10 



Table 6. The 68% and 95% confidence limits for the growth factor 
7 obtained for our paramctcrisation with DGP as the fiducial 
model compared to those obtained using only linear scales and 
compared to the constraint from using Smith ct al. 2003 to model 
the non-linear. These are marginalised over f2 m , erg, A, ai and 
Ot2 ■ 
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0.20 0.34 0.48 0.62 0.76 0.90 



0.45 0.55 0.65 



(a) ACDM fiducial model with ground-based errors 



(b) ACDM fiducial model with Euclid errors 



0.0 0.5 1.0 1.5 2.0 2.5" 




0.0 0.3 0.6 0.9 1.2 1.5 




(c) DGP fiducial model with ground-based errors 



(d) DGP fiducial model with Euclid errors 



Figure 9. Constraints on 7, cti, 02 and A from our ground-based survey and Euclid, using 0.4 redshift bins between 0.3 and 1.5, where 
we have marginalised over all Q m and erg. The light grey contours show the 68% confidence limits and the dark grey show the 95% 
confidence limits. 



and [6] show the percentage difference between the 68% con- 
straint obtained for 7 if only a linear analysis is used com- 
pared to the full non-linear analysis with the fitting formula 
is 100% for the ground-based survey and 300% for Euclid 
with a ACDM fiducial model and 70% for the ground-based 
survey and 140% for Euclid with a DGP fiducial model. 

The percentage overestimation, shown in Tables [5] and 
[6] at the 68% level, in the ability of the ground-based survey 
and Euclid to constrain 7 if only the Smith et. al. fitting 
formula is used is 10% for the ground-based survey and 40% 
for Euclid with a ACDM fiducial model and 10% for ground- 
based survey and 60% for Euclid with a DGP fiducial model. 
This demonstrates that if a full non-linear analysis is to be 



used then it is necessary to ensure that GR is obtained at 
small scales, and the extra parameters from the Hu and 
Sawicki fitting formula must also be measured. 



5 CONCLUSIONS 

In this paper we have presented weak lensing predictions for 
modified gravity models, including the non-linear regime of 
the power spectrum. 

We have shown how the power spectrum is calculated 
fo r DGP, f(R) and QCD M models, using the fitting function 
of lHu fc Sawickil i|2007al ) to explore deep into the non-linear 
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regime, while including the fact that gravities should tend 
towards GR on small scales. 

We have calculated the total shear power spectrum 
given the modified gravity power spectrum, and have shown 
that this will be measured with high signal-to-noise with fu- 
ture lensing surveys such as Euclid and DES. We have taken 
into account the cosmic covariance in addition to the noise 
due to the intrinsic shapes of galaxies. 

We have shown that there is substantial additional dis- 
criminatory power between modified gravity models which 
is now afforded to us by the inclusion of the non linear power 
regim e. We have also shown that using only the lSmith et al.l 
(2003) formula without any attempt to obtain the GR non- 
linear power spectrum on small scales leads to an overes- 
timatation in the ability of future surveys to differentiate 
between different growth histories. 

We have parameterised the dark matter power spectrum 
using the growth factor 7 and the parameters in the non- 
linear fitting function to see how well a ground-based survey 
similar to DES, and a space-based survey such as Euclid, 
will be able to put constraints on these. We have compared 
the results from this parameterisation with results obtained 
from using only linear scales and have shown the constraint 
on 7 to be much tighter in the former case. 
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